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ABSTRACT 

We report the discovery and early observations of the peculiar Type Iln supernova (SN) 2006gy 
in NGC 1260, revealing that it reached a peak visual magnitude of about —22, making it the most 
luminous supernova ever recorded. It had a very slow rise to maximum that took about 70 d and 
stayed brighter than —21 mag for about 100 d. It is not yet clear what powers the enormous luminosity 
and the total radiated energy of ~ 10 51 erg, but we argue that any known mechanism — thermal 
emission, circumstellar interaction, or 56 Ni decay — requires a very massive progenitor star. The 
circumstellar interaction hypothesis would require truly exceptional conditions around the star, which, 
in the decades before its death, must have experienced a luminous blue variable (LBV) eruption like 
the 19th century eruption of 77 Carinae. However, this scenario fails to explain the soft and unabsorbed 
X-ray emission detected by the Chandra X-ray Observatory, which constrains the progenitor's mass- 
loss rate to be three orders of magnitude too small. Alternatively, radioactive decay of 56 Ni may 
be a less objectionable hypothesis, but it would imply a large Ni mass of ~22 M Q , requiring that 
SN 2006gy was a pair-instability supernova where the star's core was obliterated (rather than forming 
a neutron star or a black hole). While this is still uncertain, SN 2006gy is the first supernova for 
which we have good reason to suspect a pair- instability explosion. Independent of the energy budget, 
a narrow Ha emission line from unshocked circumstellar gas also suggests a very massive progenitor 
star. The shell has a mass of several M of hydrogen and expansion speeds of 130-260 km s _1 , ruling 
out progenitor stars with initial masses below 40 M . Based on a number of lines of evidence, we 
eliminate the hypothesis that SN 2006gy was a "Type Ha" event — that is, a white dwarf exploding 
inside a hydrogen envelope. Instead, we propose that the progenitor may have been a very massive, 
evolved object like r\ Carinae that, contrary to expectations, failed to completely shed its hydrogen 
envelope before it died. Our interpretation of SN 2006gy implies that some of the most massive stars 
can explode prematurely during the LBV phase, preventing them from ever becoming Wolf-Rayet 
stars. SN 2006gy also suggests that some of the most massive stars can create brilliant supernovae 
instead of experiencing ignominious deaths through direct collapse to a black hole. If such a fate is 
common among the most massive stars, then observable supernovae from Population III stars in the 
early universe will be more numerous than previously believed. 

Subject headings: circumstellar matter — stars: evolution — supernovae: individual (SN 2006gy) 



1. INTRODUCTION 

Supernovae (SNe) resulting from the deaths of massive 
stars span a wide range of peak absolute visual mag- 
nitude, typically between —15 and —20.5, and usually 
reaching their peak within about 20 d. They also exhibit 
a range of spectral properties depending on the extent to 
which products of nuclear burning are exposed at their 
surface, as well as the expansion speed and the amount 
of circumstellar material. Their diversity depends on 
the star's initial mass and rate of mass loss during its 
lifetime. Current expectations are that stars born with 
initial masses above ~40 Mq, which never become red 
supergiants (RSGs; Humphreys & Davidson 1979), will 
shed their hydrogen envelopes to expose their He core 
before they die (e.g., Abbott & Conti 1987). As Wolf- 
Rayet (WR) stars, they are then expected to explode, 
producing Type Ib/c SNe (see Filippenko 1997). Based 
on observations of SN 2006gy that we discuss here, we 
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speculate that this scenario does not always apply. 

One way to prevent a star from reaching the WR phase 
before explosion would be if the star's mass-loss rate is 
insufficient to shed the hydrogen envelope before the end 
of core He burning. This is thought to be the case for 
massive stars in the early universe, because their much 
lower (or zero) metallicity should make their line-driven 
stellar winds very inefficient (Baraffe et al. 2001; Ku- 
dritzki 2002; Heger et al. 2003). Depending on the mass 
at the time of death, very massive stars in this predica- 
ment might suffer a pair-production instability explosion 
(Barkat et al. 1967; Fraley 1968; Bond, Arnett, & Carr 
1984; Heger & Woosley 2002), where the star's core is 
obliterated instead of collapsing to a black hole. 

However, there are reasons to suspect that the mass- 
loss properties of stars in the local universe may not be so 
different from these early stars. Namely, recent studies of 
line-driven winds from O-type stars and WR stars have 
shown that their winds are highly clumped, requiring 
that their mass-loss rates through line-driven winds on 
the main sequence could be an order of magnitude lower 
than previously believed (Fullerton et al. 2006; Bouret et 
al. 2005). In that case, for stars with initial masses above 
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~40 M Q that never become RSGs, the burden of mass 
loss falls to the post-main-sequence luminous blue vari- 
able (LBV) phase, when very massive stars suffer mul- 
tiple giant eruptions that shed several M Q in just a few 
years (Smith & Owocki 2006). If these LBV eruptions 
are not sufficient to remove the star's entire outer hy- 
drogen envelope fast enough, as may be the case for the 
most massive stars above 100 M Q , then the star would 
seem to explode early as an LBV producing a Type Iln 
event. Interestingly, Gal- Yam et al. (2007) find that the 
rate of Type Iln events is in broad agreement with the 
hypothesis that they are the explosions of extreme LBVs. 
The fact that giant LBV eruptions are continuum driven 
may hint that low-metallicity stars may be capable of 
shedding mass after all (Smith & Owocki 2006), which 
would affect the range of initial masses that are subject 
to the pair instability in Population III stars. Because 
stars that begin their lives above 100 M Q are so few in 
number, their end fates are poorly constrained by obser- 
vations (see Gal- Yam et al. 2007 for relevant discussion), 
and are still an open question. For these reasons, any 
potential detection of a pair-instability supernova in the 
modern universe would be of great interest to stellar as- 
trophysics. Here we explore this notion, along with oth- 
ers, as a possible explanation for the bizarre properties 
of SN 2006gy. 

SN 2006gy in the peculiar SO/Sa galaxy NGC 1260 
was discovered and confirmed by the Texas Supernova 
Search (TSS; Quimby 2006a) with the ROTSE Illb tele- 
scope (Akerlof et al. 2003) at McDonald Observatory in 
unfiltcrcd images (Quimby 2006b) taken on 2006 Sep. 
18.3 (UT dates are used throughout this paper). It was 
initially classified (Harutyunyan et al. 2006) as a SN II 
(actually SN Iln, based on the written description), but 
Prieto et al. (2006) nearly simultaneously suggested that 
the object was instead a bright active galactic nucleus 
(AGN). However, in the subsequent month, our group 
continued to follow SN 2006gy, and with additional as- 
trometric, photometric, and spectroscopic data we an- 
nounced that it did indeed appear to be a SN after all, 
and not an AGN (Foley et al. 2006). In this paper we 
present additional data and analysis of SN 2006gy, lead- 
ing us to propose that it marked the death of a very 
massive star with much of its hydrogen envelope still in- 
tact, while surrounded by a massive circumstellar neb- 
ula. In many respects, the type of progenitor we infer 
for SN 2006gy resembles the LBV star rj Carinae in our 
own Galaxy, as discussed below. 

2. OBSERVATIONS 

2.1. Imaging and Photometry 

Figure 1 shows a laser guide star (LGS) adaptive optics 
(AO) near-infrared image of SN 2006gy and the nucleus 
of its host galaxy NGC 1260, revealing a clear offset of 
the SN from the galaxy center. Images at three wave- 
bands (J, H, and K s ) were obtained on 2006 Nov. 4 
using the AO system in LGS mode (Lloyd et al. 2000; 
Max et al. 1997) on the Shane 3-m telescope at Lick Ob- 
servatory. The total integration time in each band was 
480 s, accumulated over 8 exposures. The native scale 
of the 256 x 256 pixel Rockwell PICNIC array is 0"076 
pixel -1 (Perrin 2007). Mosaiced images have a scale of 
0V04 pixel -1 . The SN itself was bright enough to use as 
a "tip-tilt" star for the LGS system. The effective res- 
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Fig. 1. — Laser guide star adaptive optics image of SN 2006gy 
and the nucleus of NGC 1260, showing a clear offset of the SN from 
the galaxy center. Blue is J band (1.25 /im), green is H band (1.65 
/im), and red is K s band (2.2 fim). 

olution (full width at half-maximum intensity; FWHM) 
is 0'.'2 in the H band. The measured offset of the SN 
from the centroid of the galactic nucleus is 0'.'941 west, 
0'.'363 north, with a la uncertainty of O'.'Ol in each direc- 
tion; this confirms and improves the earlier offset mea- 
surement (Foley et al. 2006) of 0'.'880 west, 0'.'140 north, 
±0'/08. SN 2006gy is therefore located about 350 pc from 
the galaxy's center (at its assumed distance of ~73 Mpc), 
confirming that it is not an AGN. 4 

Figure 2 shows the i?-band light curve of SN 2006gy 
obtained by our group using the Katzman Automatic 
Imaging Telescope (KAIT; Filippenko 2003) at Lick Ob- 
servatory, compared to a sample of several other rep- 
resentative SN light curves. The unfiltered KAIT im- 
ages for SN 2006gy were used to derive an U-band light 
curve. As demonstrated by Riess et al. (1999) and Li 
et al. (2003), the best match to broad-band filters for 
the KAIT unfiltered data is the R band. Each image is 
aligned to a deep pre-SN image, and the contamination 
of the host-galaxy emission is carefully removed. The 
net flux for the SN is then compared to 19 bright stars 
using calibrations from the USNO Bl catalog. We list 
the KAIT apparent R magnitudes of SN 2006gy in Table 
1. To put the flux of SN 2006gy on an absolute magni- 
tude scale, we adopt a distance to the host galaxy NGC 
1260 of 73.1 Mpc, using H = 72 km s -1 Mpc" 1 and 
using a recession velocity for the central cluster galaxy 
of 5361 km s . We also assume a Galactic reddening of 
Ar = 0.43 mag (Schlegel et al. 1998) and a host-galaxy 
reddening of An — 1.25 ± 0.25 mag (see §2.2 and Fig. 
3). In Figure 2 we plot days since explosion instead of 
days since discovery. Our first measurement with KAIT 
was a non-detection made on 2006 Aug. 26, which was 
23 d before the discovery of SN 2006gy. Judging from 
the slowly rising curve, we estimate that the explosion 
date was roughly 6 d before the KAIT non-detection. 

2.2. Lick and Keck Spectroscopy 

Figure [3] shows two visual-wavelength spectra of 
SN 2006gy obtained on 2006 Sep. 25.5 and 2006 Oct. 

4 Ironically, NGC 1260 may contain a faint AGN after all, al- 
though SN 2006gy is a real SN explosion. Later in this paper we 
also present an X-ray image of SN 2006gy which shows two sources, 
one being the SN and the other the nucleus of NGC 1260. 
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KAIT Photometry of SN 2006gy 
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30.4 using the Kast double spectrograph (Miller & Stone 
1993) mounted on the Lick Observatory 3-m Shane tele- 
scope. The long slit of width 2" was aligned along the 
parallactic angle (Filippcnko 1982). The data were re- 
duced using standard techniques as described by Foley 
et al. (2003) and references therein. The spectra were 
corrected for atmospheric extinction (Bessell 1999; Math- 
eson et al. 2000) and then flux calibrated using standard 
stars observed at an airmass similar to that of the SN. 

The closest match to SN 2006gy in our spectral 
database is SN 2006tf, taken on 2007 Jan. 13, as shown 
in Figure [3] 5 The red continuum shape of SN 2006gy 
is unusual for SNe Iln, which are typically much bluer 
(Schlegel 1990), so we have plotted the SN 2006gy spec- 
trum after removal of various amounts of reddening for 
comparison. Although a direct comparison to SN 2006tf 

J SN 2006tf was discovered in the course of the Texas Supernova 
Search on 2006 Dec. 12 UT (Quimby et al. 2007). With a discovery 
magnitude of 16.7 and a redshift z = 0.074, the SN has an absolute 
magnitude of —20.7, which is extremely luminous but still less so 
than SN 2006gy. Our follow-up photometry also suggests that 
SN 2006tf exhibits a light-curve shape similar to that of SN 2006gy. 




Days since explosion 



Fig. 2. — Comparison of the absolute ij-band light curve of 
SN 2006gy with those of other SNe. We plot days since explosion, 
which we judge to be ~29 d prior to the discovery of SN 2006gy. SN 
1998dh is a typical SN la, and the data are from our unpublished 
photometric database, with a typical absolute magnitude of Mr = 
— 19.5 mag assumed. SN 1999em is a typical Type II (Leonard et 
al. 2002), SN 19941 is a well-observed SN Ic (Richmond et al. 1996), 
and SN 1998bw is a peculiar SN Ic (Galama et al. 1998). SN 1987A 
is a peculiar SN II, with a broad light curve but a low luminosity 
(from Hamuy et al. 1990). SN 1994W is a SN Iln that is powered 
by strong interaction with its circumstellar material (Sollerman et 
al. 1998). We also plot two unusual SNe that are relevant to the 
discussion of SN 2006gy: SN 2002ic (Hamuy et al. 2003) and SN 
2005gj (Aldering et al. 2006). 

is complicated by the temporal evolution, the early (day 
36) spectrum of SN 2006gy seems most consistent with 
Ar = 1.5 mag, while the later (day 71) spectrum is 
more consistent with Ar = 1.0 mag (the spectra of 
SN 2006gy were already corrected for Galactic extinc- 
tion of Ar = 0.43 mag, as noted earlier). Comparison 
with other SNe Iln at similar phases (not shown) also 
suggests values of the host-galaxy value of Ar — 1.0 to 
1.5 mag. We therefore adopt Ar = 1.25 ± 0.25 mag for 
SN 2006gy. The extinction could be higher if SN 2006tf 
has its own significant reddening, although it appears to 
have very weak Na I D absorption. The strong Na I D 
absorption in the spectrum of SN 2006gy may suggest 
higher reddening than we have assumed here, so our es- 
timates of luminosity for SN 2006gy are conservative. 

Figure [4] shows the day 36 Lick spectrum from Fig- 
ure [H and also a spectrum with a smaller wavelength 
range and higher spectral resolution of R « 4500 taken 
near maximum light on day 96. The latter spectrum 
was obtained on 2006 Nov. 24.51 using the DEIMOS 
spectrograph (Faber et al. 2003) on the Keck II tele- 
scope. Using a customized version of the DEEP data 
reduction pipeline, we obtained sky-subtracted, recti- 
fied two-dimensional images, and wavelengths were cali- 
brated with respect to an internal calibration lamp (Foley 
et al. 2007). We checked carefully to make sure that the 
sky-subtraction procedure did not artificially introduce 
narrow absorption components; this is implausible based 
on the final results anyway, since H and He I lines show 
similar blueshifted absorption profiles. We corrected for 
telluric absorption (Matheson et al. 2000) by comparison 
with the standard star BD+284211. 

Figure \E\ shows the Ha profile of SN 2006gy near max- 
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Fig. 3. — Lick Observatory spectra of SN 2006gy at two different epochs, corrected for a range of assumed host-galaxy reddening 
corresponding to the values of An listed at right (Cardelli et al. 1989). This extinction is in addition to Galactic extinction of An = 0.43 
mag. These are compared to the day 32 spectrum of the Type Iln SN 2006tf (black) from our database, which is a SN with a spectrum 
similar to that of SN 2006gy, but seems to show little reddening. We adopt An = 1.25 ± 0.25 mag for SN 2006gy; see text. 
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Fig. 4. — Dereddened visual-wavelength spectra of SN 2006gy at t = 36 d and 96 d after explosion, obtained at Lick Observatory and 
with the Keck II telescope, respectively. Several narrow absorption lines in our high-resolution Keck spectrum have been marked, but 
there are some remaining unidentified lines. Also plotted is a spectrum of the Type la SN 1991T at t = 35 d (Filippenko et al. 1992) for 
comparison with our day 36 spectrum of SN 2006gy; there is essentially no similarity between the two spectra. 



imum light from a portion of the same Keck spectrum in 
Fig. |4l with the flux normalized to the underlying contin- 
uum level, and the velocity scale chosen with the narrow 
Ha emission feature at v — km s _1 . The Ha profile in 
Figure \5\ reveals several different characteristic velocities 
relevant to interpretations of SN 2006gy. First, the very 
narrow emission component (FWHM » fOO km s _1 ) has 
an associated P Cygni absorption feature that indicates 
outflow speeds of 130 km s" 1 (the trough) to 260 km s _1 
(the blue edge) in the unshocked circumstellar gas. In 
addition to Ha, several lines identified in Figs |4] and [5] 
also have narrow absorption features. 

A broad Ha emission component has an apparent 
FWHM « 2400 km s" 1 that is similar to H/3 at early 
times (Harutyunyan et al. 2006). The true unabsorbed 
FWHM of this broad Ha component is larger because of 
the broad blueshifted absorption. Extended faint wings 



out to ±6,000 km s _1 may be caused either by electron 
scattering or by the fastest SN ejecta. 

The blue edge of the broad, blueshifted Ha absorption 
in Figure [5] indicates an outflow speed of 4,000 km s — 1 , 
where the emission jumps back up just to the level that 
would be expected for a symmetric profile. This jump is 
readily apparent when we take the redshifted side of the 
broad emission profile and reflect it to the blue side, to 
simulate what a symmetric profile would look like (Fig. 
[5]) . Because this absorption traces the speed of the dom- 
inant absorbing material along the line of sight at this 
epoch, we take this speed of 4,000 km s _1 to represent 
dense material swept up by the SN blast wave in the cir- 
cumstellar material (CSM) interaction hypothesis, which 
should closely trace the speed of the blast wave itself. 

The broad-line profile differs from the smooth broad 
parts of Ha profiles normally seen in SNe Iln (e.g., 
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Chugai et al. 2004). The blueshifted absorption trough 
flattens out and does not descend below the underlying 
continuum level. This may hint that the continuum lu- 
minosity and Ha emission/absorption have different ori- 
gins, and provides important clues to the shell optical 
depth and CSM density. For example, the blueshifted 
absorption may arise in shocked CSM gas, whereas the 
continuum luminosity may originate in the SN ejecta. 
Asymmetric geometry in the CSM obviously may be rel- 
evant. These details have some bearing on the hypothe- 
ses for the power sources discussed in §3.3 and §3.4. In 
any case, this broad, blueshifted Ha absorption probably 
shares an origin with the broad, blueshifted absorption 
features for other lines identified in Figure [U In light 
of possible geometric complexities, we defer a detailed 
discussion of the line profiles to a later paper. 

2.3. X-ray Observations, Data Reduction, and Analysis 

The Chandra X-ray Observatory began observing the 
location of SN 2006gy on 2006 Nov 14.86 using Direc- 
tor's Discretionary Time. The observation lasted 29.743 
ks, and the data were taken with the Advanced CCD 
Imaging Spectrometer using an integration time of 3.2 s 
per frame. The telescope aimpoint was on the back-side 
illuminated S3 chip, and the data were telemetered to 
the ground in "very faint" mode. 

Data reduction was performed using the CIAO 3.4 soft- 
ware provided by the Chandra X-ray Center 6 . The data 
were reprocessed using the CALDB 3.3.0 set of calibra- 
tion files (gain maps, quantum efficiency, quantum effi- 
ciency uniformity, effective area) including a new bad- 
pixel list made with the acis_run_tiotpix tool. The re- 
processing was done without pixel randomization that is 
added during standard processing. This omission slightly 
improves the point-spread function (PSF). The data were 
filtered using the standard ASCA grades (0, 2, 3, 4, and 
6) excluding both bad pixels and software-flagged cosmic- 
ray events. A search was done for strong background 
flaring, but none was found. 

Absolute Chandra astrometry is typically good to 0'.'5, 
and we sought to tie the Chandra frame to the KAIT 
image to obtain a reliable identification of the nucleus 
of NGC 1260 and the SN in the Chandra data. Sev- 
eral Chandra point sources were found using the CIAO 
wavdetect tool, and their positions were refined using 
ACIS Extract version 3.107 (Broos et al. 2002). Three of 
these sources had KAIT counterparts, although one had 
a somewhat poorly determined Chandra position due to 
its location ~3' off-axis (the Chandra PSF degrades as 
a function of off- axis angle). Using all three sources, we 
obtained an astrometric correction to the Chandra data 
of 0'.'329 in right ascension (a) and 0'.'089 in declination 
(S). Using the two best counterparts, we obtained shifts 
of Aa = 0'.'472 and AS = 0'.'104. We use this latter shift 
for the rest of our analysis. 

Figure [6] shows a 0.5-2 keV image of the Chandra data 
after this shift; arrows indicate the KAIT positions of the 
SN (red) and galaxy nucleus (blue). In addition to the 
raw image, Figure [6] shows a Gaussian-smoothed image 
and a maximum-likelihood reconstruction of the data, as 
well as an image of the Chandra PSF on the same spatial 
scale. The maximum-likelihood reconstruction was made 



by ACIS Extract using the max_likelihood procedure 
available in the IDL Astronomy User's Library 7 ; we went 
through 200 iterations of the algorithm, using the PSF 
shown in the figure. The PSF was constructed by ACIS 
Extract through use of the CIAO tool mkpsf based on 
the off-axis location of the source and at an energy of 
1.49 keV (the Chandra PSF is also a function of energy). 
As can be seen, there is excellent agreement between the 
locations of the reconstructed sources and the locations 
of the SN and host-galaxy nucleus. This argues strongly 
that we have, in fact, detected SN 2006gy and spatially 
resolved it from the nucleus of NGC 1260. 

We measured counts in the full 0.5-8 keV bandpass 
from the position of the SN using a small extraction re- 
gion to minimize contamination from the galaxy nucleus. 
The extraction region has a radius of ~0 / /4, correspond- 
ing to about 40% of the PSF. Response files were con- 
structed with the CIAO tools, and ACIS Extract cor- 
rected them for the non-standard extraction region. The 
background region is a source-free annulus centered on 
the position of the SN with inner and outer radii of 6" 
and 14", respectively. Based on the 241 counts detected 
in this region, we expect only 0.24 background counts in 
our extraction region. In the restricted energy range of 
0.5-2 keV (used for the rest of this paper), we expect 
only 0.08 background counts in our extraction region. 

Four counts were detected in our extraction region, 
which precludes a detailed spectral analysis. However, 
the counts were all detected below 2 keV, giving some 
indication of the spectral shape. We assume a thermal 
plasma spectrum (Raymond-Smith) with kT — 1 keV 
to estimate the luminosity. Such thermal spectra have 
successfully fit the X-ray spectra of SNe, and tempera- 
tures much higher than this would result in significant 
emission detectable by Chandra (which was not seen). 
Based on an assumed reddening toward SN 2006gy of 
E(B — V) = 0.74 mag, we assume an X-ray absorbing 
column of n# = 4.1 x 10 21 cm" 2 (Predehl & Schmidt 
1995). Such an absorbed thermal plasma observed by 
Chandra would result in a ratio of 0.5-2 keV to 2-8 keV 
counts of ~10:1, in accordance with observations. We fit 
this model to the observed 0.5-8 keV spectrum in Sherpa 
(Freeman et al. 2001) using the statistic of Cash (1979). 
The only free parameter is the overall normalization of 
the model. From the best fit we find an unabsorbed X- 
ray luminosity (0.5-2 keV) of 1.65 x 10 39 erg s _1 . 

3. THE DEATH OF A VERY MASSIVE STAR WITH ITS 
HYDROGEN ENVELOPE INTACT 

3.1. The Energy Budget and a High-Mass Progenitor 

SN 2006gy has quickly distinguished itself as unique 
from other SNe in two important ways. First, after cor- 
recting for distance and extinction, it is the most lumi- 
nous SN ever seen, and second, it has exhibited a re- 
markably slow rise to its peak luminosity, and has stayed 
bright for a very extended time. SN 2006gy has peaked 
and is now on a slow decline, but even after 200 d it is 
still as luminous as the peak of a typical SN la. 

SN 2006gy was classified as a SN Iln with narrow hy- 
drogen lines in its spectrum at early times (Harutyunyan 
et al. 2006), although the spectrum has notable differ- 
ences compared with prototypes of this class. It dramat- 
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ically violates the expectation that SNe II are generally 
less luminous than SNe la (Fig. 2 includes a fairly typical 
Type II SN 1999em), and that SNe Iln usually take only 
~20 d to reach their peak (Li et al. 2003). SN 2006gy, 
by contrast, took ~70 d to gradually climb to its peak. 
For about 100 d it was more luminous than Mr = — 21 
mag, brighter than any other SN known to date. 

Simply put, for a supernova to be extremely luminous 
and to remain that way for such an extended time is truly 
spectacular. Integrating the light curve in Figure 2 and 
assuming zero bolomctric correction, we calculate a total 
radiated energy of £? rad = (1.2 ± 0.2) x 10 51 erg. This 
requires either very efficient conversion of blast-wave ki- 
netic energy into light, or some alternative energy source. 
One or a combination of the three following traditional 
mechanisms may power the visual light: (1) H recom- 
bination/thermal radiation of the supernova ejecta, (2) 
interaction of the supernova blast wave with the CSM, or 
(3) energy from radioactive decay of 56 Ni. Continued ob- 
servations and probably extensive theoretical work will 
be needed to choose decisively between these options, 
but here we argue that regardless of which of these three 
mechanisms is responsible, the extreme energy budget of 
SN 2006gy requires that its very massive progenitor star 
retained its H envelope until it exploded. 

The first option of thermal emission from the H- 
recombination front in the supernova debris would re- 
quire a huge ejected mass of order 100 M© or more, based 
simply on the total radiated energy. A heavy H enve- 
lope might help explain the unusually slow speed of only 
about 4000 km s _1 indicated by the Ha line (Fig. [5]), and 
might provide a natural explanation for the long dura- 
tion and rise time of the SN because of time needed for 
energy to diffuse out of the massive envelope. Whether 
or not the SN could actually radiate efficiently enough 
to produce the observed luminosity with this mecha- 
nism remains to be proven and should be investigated 
with detailed calculations. For example, at the tempera- 
ture of the photosphere defined by the H-recombination 
front (typically 5000 K to 8000 K), the luminosity of 
SN 2006gy requires an emitting radius larger than what 
we might expect from its observed expansion speed of 
4000-4500 km s" 1 and age. Instead of 70 d, the ob- 
served peak luminosity would seem to require an age of 
200-380 d since explosion (assuming linear motion), or 
rapid deceleration at early times. Such rapid decelera- 
tion at early times cannot be ruled out by our data. 

The second option of powering the visible light entirely 
with CSM interaction is problematic, but is difficult to 
rule out conclusively. ^From the relatively weak soft X- 
ray flux of SN 2006gy detected by Chandra, we derive 
an upper limit to the progenitor star's mass-loss rate of 
- 5 x 10~ 4 M yr" 1 (see §3.2). We find that this falls 
short of the circumstellar density that would be needed 
to power the visual light curve of SN 2006gy by three 
orders of magnitude (§3.3). In order to explain the high 
luminosity in whole or in part by CSM interaction, one 
would therefore need to assume that the X-ray emission is 
severely quenched and that the Chandra detection is er- 
roneous — but this is difficult to accept, since we clearly 
detect soft (unabsorbed, not hard) X-ray emission from 
the position of the SN (Fig. [6|. Even if it were true, 
though, a closer look at the demands placed on the cir- 
cumstellar density make it difficult to explain with any- 
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Fig. 5.— The Keck/DEIMOS spectrum of the Ha line seen in SN 
2006gy, with the flux normalized to the underlying continuum. The 
upper-right inset shows a closer view of the narrow P Cygni line 
profile that we believe to be associated with dense unshocked CSM. 
The blueshifted narrow absorption trough has a minimum at about 
— 130 km s — 1 , reaching —260 km s _1 at its blue edge. The other 
narrow absorption lines labeled as "Fe II" are Fe II AA6418, 6433, 
6456, 6517. The dashed line labeled "symmetric" is the red side 
of the broad Ha line reflected to blueshifted velocities, showing 
what the line shape would be if it were symmetric. Comparing 
this to the observed Ha profile, we see significant blueshifted Ha 
absorption from km s _1 out to a sharp blue edge at about —4000 
km s , which we take to be the dominant speed of the SN blast 
wave. At that point, the blueshifted emission recovers to the level 
expected for a symmetric profile, and then gradually declines to 
the continuum level at about —6000 km s - , just as on the red 
side of the line (which overlaps with He I A6680). 

thing other than a massive star that coincidentally had 
an LBV outburst just before the supernova explosion. 

Finally, the third option, radioactive decay of 56 Ni, is 
perhaps the least problematic, as we will discuss further 
in §3.4. The main point of interest is that if this mech- 
anism powers the visual light, then the high luminosity 
of SN 2006gy requires a very large Ni mass that cannot 
arise from a normal core-collapse SN. Instead, the large 
mass involved would require that SN 2006gy was a pair 
instability supernova in which the star's core was oblit- 
erated. If true, SN 2006gy would be the first observed 
example of a pair-instability supernova. This mechanism 
also has some potential difficulties, but they are more 
along the lines of uncharted theoretical territory, rather 
than fundamental physical or observational constraints. 
Therefore, SN 2006gy provides fertile ground for impor- 
tant theoretical work in this area. 

3.2. Limits to the Progenitor's Mass-Loss Rate from 
X-ray Data 

If we interpret the X-ray emission detected by Chandra 
as the result of interaction of the outgoing shock with 
circumstellar material (CSM interaction), we can place 
an upper limit on the mass-loss rate of the progenitor 
star. This interaction has been explored in detail (e.g., 
Fransson, Lundqvist, & Chevalier 1996). The softness of 
the X-ray emission points toward a reverse-shock origin, 
and we use the adiabatic case. A useful form of their eq. 
(3.10) is found in Pooley et al. (2002): 
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= 2 x 10 35 C(n - 3)(n - 4) 2 r^ - 24 ^ - 116 /^ 

Kfe) 2 ^ 1 ^)" 1 ergs^keV" 1 , (1) 

where £ = 0.86 for solar abundances, n is the index of 
the eject a- density profile [psn oc t~ 3 (r/t)~ n ], T 8 is the 
temperature in units of 10 8 K, M_6 is the progenitor's 
steady-state mass-loss rate in units of 10~ 6 Mq yr^ 1 , 
V w \ is its wind speed in units of 10 km s" 1 , V S 4 is the 
shock velocity in units of 10 4 km s _1 , and t is the time 
since explosion. 

The value of n appropriate for SN 2006gy is uncertain, 
but typical values for core-collapse SNe are in the range 
7-12. We assume a temperature of 1 keV, for which 
Tg = 0.116. From Figure [5] we take the wind speed to 
be ~200 km s _1 (V w i — 20), and the shock velocity to 
be 4500 km s _1 (V S 4 = 0.45). The Chandra observation 
took place 87 d after the explosion. 

This implies a mass-loss rate for the progenitor of 
1.4 x 10 -4 M© yr" 1 assuming a steady mass- loss rate 
of the progenitor in the decades before explosion, and 
adopting a SN ejecta density profile with n = 12. For 
a profile with n = 7, the mass-loss estimate rises to 
5.4 x 10 -4 M Q yr" 1 . This range of mass-loss rates is 
in good agreement with observed values in luminous fi- 
nch WN stars (e.g., Hamann et al. 2006) or quiescent 
non-outburst LBVs (Smith, Vink, & de Koter 2004). 
As we discuss below, however, this range of mass-loss 
rates falls short of that needed to power the luminos- 
ity of SN 2006gy with CSM interaction by three orders 
of magnitude. This is a serious obstacle to any such 
model, which must now account for why we observe a rel- 
atively weak and soft (i.e., unabsorbed) X-ray flux from 
SN 2006gy. A likely explanation is that CSM interaction 
is important in creating the observed soft X-rays and in 
causing the emission-line spectrum of SN 2006gy (espe- 
cially the broad Ha emission), but that something else 
drives its visual-wavelength continuum luminosity. Be- 
low, we consider the CSM interaction hypothesis (§3.3) 
as a power source for SN 2006gy aside from the difficulty 
posed by X-rays, as well as an alternative energy source 
for its radiated luminosity (§3.4). 

3.3. A Closer Look at Circumstellar Interaction 

Ofek et al. (2007) suggested CSM interaction as a 
means to power the visual light of SN 2006gy, but here 
we wish to make a clear distinction between two differ- 
ent scenarios. The first is where the blast wave from 
a SN la interacts with dense CSM from a companion 
star that provides the hydrogen in the spectrum (the so- 
called "Type Ha" scenari o; e.g., Deng et al . 2004), as 
suggested in version 1 (in astro-ph/06 12408) of the re- 
cent study by Ofek et al. (2007). This interpretation 
had also been suggested previously for the bright SNe 
2002ic and 2005gj (Hamuy et al. 2003; Deng et al. 2004; 
Aldering et al. 2006). Note, however, that Benetti et al. 
(2006) have instead argued in favor of a core-collapse ori- 
gin for SN 2002ic, so the true nature of these events is still 
controversial. The second type of scenario would be a 
blast wave from a core-collapse or pair-instability super- 
nova from a massive star interacting with its own ejecta, 



analogous to the interpretation of the SN Iln 1994W by 
Chugai et al. (2004). 

We argue here that the first scenario (SN Ha) is un- 
tenable for SN 2006gy for a number of reasons. Based 
in part on a preprint of our work presented here (version 
1 of |astro-ph/0612617| ), Ofek et al. (2007) revised their 
original Type Ha interpretation of SN 2006gy to include 
the possibility that it could have been a massive star as 
we originally proposed. The second SN 1994W-like sce- 
nario, on the other hand, is almost certainly relevant to 
SN 2006gy, but based on the weak X-ray emission we 
probably require a different source for the bulk of the 
radiated luminosity. If, for the sake of argument, we 
demand that CSM interaction powers the luminosity, we 
find that the extraordinary energy demands of SN 2006gy 
point to a circumstellar environment that is only likely 
to be produced by a very massive star that suffered a 
rare outburst immediately prior to the SN. In the case 
of SN 2006gy, the luminosity and total energy need to 
be scaled up by a factor of 40 or more from those for 
SN 1994W. 

In order to power the luminosity of SN 2006gy with 
CSM interaction, the environment created by the progen- 
itor star must be extraordinarily dense. Ofek et al. (2007) 
originally (version 1 of astro-ph/0612408) estimated that 
to achieve the luminosity of SN 2006gy with a shock 
plowing into CSM, the progenitor star (or its compan- 
ion star in a close binary system) needed to have a wind 
with an average mass-loss rate of ^10~ 2 M Q yr _1 in the 
decades before explosion. However, this estimate scales 
with the adopted wind speed V w and inversely with the 
shock speed V s , which Ofek et al. originally took to be 
V w = 10 km s _1 and V s — 10 4 km s _1 . Instead, though, 
we observe a much faster speed of V w rj 200 km s _1 in 
the circumstellar environment indicated by the narrow P 
Cygni component in our spectra (Fig. [5j see §3.5), rais- 
ing this necessary mass- loss rate to ~ 0.2 M Q yr _1 to 
achieve the same circumstellar density (Ofek et al. 2007 
note this in version 2 of their paper, based on velocities 
in our Fig. [5j. We also see a slower speed for the SN 
shock of only V s w 4000 km s" 1 (Fig. || instead of 10 4 
km s _1 (Ofek et al. 2007), raising the required progeni- 
tor mass-loss rate even further to about 0.5 M Q yr -1 . 8 
Thus, if CSM interaction is to power the visual light of 
SN 2006gy, the progenitor was probably an extremely 
massive star. Recall, however, that this required value of 
0.5 M Q yr -1 is 1000 times above the highest likely value 
indicated by X-ray emission, making it problematic (see 
§3.2). Let us put this last issue aside for the time being, 
assuming that the X-rays are somehow absorbed with- 
out hardening the spectrum, so that we can consider the 
implications of the CSM interaction hypothesis. 

The expansion speed indicated by the Ha line (Fig. [5]) 
is critical for addressing the extent to which interaction 
with CSM may power the observed radiation, because the 
FWHM«2400 km s _1 of the main intermediate- width 
emission component in Figure [5] has changed little from 
the initial value of FWHM«2500 km s _1 seen in the H/3 
emission feature only a few days after discovery (Haru- 

8 One might suspect that even this value may underestimate 
what is required to power SN2006gy. In a more detailed analysis 
of SN 1994W, Chugai et al. (2004) required a similar progenitor 
mass-loss rate of 0.2 Mq yr" 1 for a short time preceding the SN, yet 
SN1994W was more than 10 times less luminous than SN 2006gy. 
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Fig. 6. — Soft-band (0.5—2 keV) Chandra images of NGC 1260. Panel a shows the raw Chandra data (after our astrometric correction) 
with red and blue arrows indicating the KAIT positions of the SN and galaxy nucleus, respectively. Panel b is a Gaussian-smoothed version 
of this image, in which the sources are more clearly apparent. Panel c is a maximum-likelihood reconstruction of the 0.5-2 keV image (see 
text for details). Panel d shows the Chandra PSF at the location of the galaxy on the same spatial scale as the other panels. 



tyunyan et al. 2006). (Recall that if the SN is pow- 
ered by CSM interaction, then the observed expansion 
speed traces the blast-wave speed, and not the decrease 
in speed expected as the H recombination front progress 
deeper into the SN ejecta.) If the expanding blast wave 
has only slowed by about 10% in the first few months, 
conservation of momentum dictates that the mass of 
swept-up material is only about 10% of the ejected mass. 
Since at least a few M of material needs to be swept up 
to power the luminosity of SN 2006gy 9 , the mass of the 
SN ejecta then needs to be at least 25 Mq. This clearly 
rules out a Type la event. Another way to approach 
the problem is that if the ejecta only slow by 10% after 
discovery, then only ^20% of the initial kinetic energy 
can be converted into radiation during that time. The 
huge radiated energy of SN 2006gy would then require 
a SN with > 5 x 10 51 erg, again too great a demand 
for a SN la, even in a double-degenerate scenario or a 
super-Chandrasekhar-mass white dwarf. 10 In short, one 
cannot extract enough energy from the shock to power 
the light curve without slowing down the shock, unless 
the initial mass and kinetic energy of the SN ejecta are 
high. 

Even if we somehow allow for very efficient conversion 
of all the 10 51 erg of blast-wave kinetic energy into ra- 

9 This comes from the required progenitor mass-loss rate, the 
duration of the SN at the time the spectrum in Fig. 5 was taken (t fs 
96 d), and the relative speed of the blast wave and circumstellar 
material: M = M X t(Vs/V w ), which gives about 2.5 Mq. 

10 Invoking the hypothesis that the CSM interaction occurred 
before the first observation, allowing the observed SN expansion 
speed to remain constant, does not help because it cannot account 
for how the light curve is powered continually for more than 100 d 
after that interaction (the ejecta cool quickly). 



diation, we must ask: What type of progenitor star is 
likely to have had such a stupendous mass-loss rate? A 
rate of 0.5 M Q yr _1 would be unheard of for a low-mass 
(2-8 M Q ) asymptotic giant branch (AGB) star, which 
is the most likely type of star to expect in the SN Ha 
scenario, for which observed mass-loss rates are four to 
five orders of magnitude lower (de Jager et al. 1988). 
Even the most extreme OH/IR stars have rates below 
10~ 4 M yr" 1 (Netzer & Knapp 1987), while the high- 
est rates during the final and brief protoplanetary neb- 
ula phase reach only (1-2) x 10~ 4 M© yr -1 (Bujarrabal 
et al. 2001). In fact, it is also more than 4 orders of 
magnitude larger than the Eddington accretion rate for 
a white dwarf, which would be relevant in a common- 
envelope scenario. Even massive stars in their normal 
(i.e., non-eruptive) states do not come close to this rate. 

The only type of star known to have a mass-loss rate 
higher than 0.1 M Q yr _1 would be an LBV during a 
giant eruption (Smith & Owocki 2006). Those events 
typically last about a decade or less (Van Dyk 2005), 
which would be of the right order (t = 150 d x Vs /V w ) to 
account for the required circumstellar environment of SN 
2006gy. The outbursts are impulsive, so the large masses 
in their nebulae (Smith & Owocki 2006) averaged over 
the durations of the visible eruptions yield these mass- 
loss rates. If it were the case that the pre-SN mass-loss 
event before SN 2006gy was of such short duration, then 
we would predict the luminosity of SN 2006gy to soon 
plummet rapidly to the late-time luminosity of a normal 
SN II. If such a drop is not observed, it will strengthen the 
case for the pair-instability hypothesis discussed next in 
§3.4. Such a sudden drop was clearly seen in SN 1994W 
at roughly day 110 (Chugai et al. 2004). 
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This interpretation, though, forces us back once again 
to the hypothesis that the progenitor was an extremely 
massive star, since only the most powerful LBV outbursts 
from the most massive stars with initial masses above 
~100 M & are known to have such high mass-loss rates. 
Coincidentally, the mass-loss rate of rj Carinae during 
its phenomenal 1843 eruption was about 0.5 M Q yr _1 
if averaged over 20 years (Smith et al. 2003). Another 
such extreme case is SN 1961V in NGC 1058 (Goodrich 
et al. 1989; Filippenko et al. 1995; Van Dyk et al. 2002), 
which is thought to have had an initial mass well above 
100 Mq. To expect such an extraordinary feat from a 
low-mass or intermediate-mass star is unreasonable even 
in the most imaginative circumstances. 

Further difficulties for the SN Ha scenario — and even 
for moderately massive progenitors — arise if we con- 
sider geometry. If one attempts to account for the un- 
usually dense circumstellar environment by invoking a 
high mass-loss rate tidal stripping "event" in a close bi- 
nary or common envelope/merger 11 , for example, then 
this would almost certainly distributed material in a flat- 
tened disk as mass is shed from the system through the 
outer Lagrangian point (e.g., Taam & Ricker 2006). In 
that case, however, even with 100% efficiency in the lo- 
cal conversion of kinetic energy into radiation, the global 
fraction of energy available is only that of the solid angle 
that can be intercepted by the disk — which will proba- 
bly be less than 10%. 

Altogether, then, there are several clear reasons why 
the Type Ha scenario originally advocated by Ofek et 
al. (2007, version 1 of |astro-ph/ 0612408) fails to power 
SN 2006gy through CSM interaction. It is perhaps not 
surprising, then, that the visual spectrum of SN 2006gy 
does not resemble a SN la or the other SN Ha candi- 
dates. Hamuy et al. (2003) argued that SN 2002ic was 
a variant of the SNe la phenomenon on the basis of the 
similarity of its spectral evolution to that of a diluted 
version of SN 1991T (Filippenko et al. 1992). While the 
continuum of the earliest spectrum of SN 2005gj was rel- 
atively featureless, it too developed the prominent broad 
iron lines typical of a SN la by two months after explo- 
sion (Aldering et al. 2006). Our earliest spectrum of SN 
2006gy is plotted in Figure H along with SN 1991T at 
a similar epoch relative to explosion. The only strong 
spectral feature in the SN 2006gy spectrum is Ha. The 
weaker features that are present do not match those of 
SN 1991T. In particular, the deep minima in the SN 
1991T spectrum near 5700 and 6200 A are lacking in 
SN 2006gy. At no later epoch did SN la features become 
visible in SN 2006gy, as can be seen in the day 71 and 
96 spectra plotted in Figures 3 and 4. We therefore have 
no compelling reason to believe that an exploding white 
dwarf was present in this event. 

We find that conversion of the blast-wave kinetic en- 
ergy into radiated luminosity might potentially power 
SN 2006gy, as has been proposed for SN 1994W (Chugai 
et al. 2004) , but only if the swept-up environment is con- 
sistent with extreme environments observed around the 
most massive evolved stars known, such as rj Carinae. 
This agrees with the conclusions in §3.5, where the prop- 
erties of the circumstellar nebula independently rule out 

11 Ignore for the moment that this hypothetical event needs to 
be synchronized with the supernova. 



progenitor stars with initial masses below 40 M Q ; initial 
masses above 60-80 M Q are favored. 

This last conclusion about the progenitor and its en- 
vironment should not be taken lightly. It requires that 
an extremely rare event analogous to the 19th-century 
eruption of rj Carinae occurred a decade or so before the 
SN explosion. Why would these two events be synchro- 
nized? We are left with a choice: Either this is such 
an unlikely event that the underlying power source for 
SN 2006gy must be some other mechanism and CSM in- 
teraction only contributes a fraction of the radiated en- 
ergy (see §3.4), or instead, it is an indication that giant 
LBV eruptions may be a sign of things to come — i.e., 
an "early warning sign" of an impending SN. The second 
possibility would be astounding if true, and SN 2006gy 
may not be alone in this regard. SN 1994W (Chugai 
et al. 2004; Sollerman et al. 1998), SN 2001em (Chugai 
& Chevalier 2006), and SN 2006jc (Foley et al. 2007) 
all show signs of dense environments that were probably 
produced by a giant mass-loss event just before the SN. 
Smith & Owocki (2006) have noted several other cases 
as well. SN 2006jc, in particular, was even observed as a 
"supernova imposter" 2 years before the final explosion 
(Nakano et al. 2006; Foley et al. 2007; Pastorello et al. 
2007). Furthermore, such an outburst preceding the SN 
event may have some theoretical expectation (e.g., Heger 
& Woosley et al. 2002). This may be a profound clue to 
the fates of the most massive stars. 

In any case, it is a marked difficulty for the CSM 
interaction hypothesis in general that — in addition 
to the softness and faintness of the detected X-rays 
noted above — the light curve, spectrum, and multi- 
wavelength properties of SN 2006gy differ from those 
of other SNe Iln powered by CSM interaction, such as 
SNe 1988Z (Filippenko 1991; Stathakis & Sadler 1991; 
Turatto et al. 1993), 1995N (Fox et al. 2000; Fransson et 
al. 2002), and 1998S (Leonard et al. 2000; Pooley et al. 
2002). SN 1988Z was bright in X-ray and radio emission 
(Schlegel & Petre 2006; Van Dyk et al. 1993; Williams 
et al. 2002), unlike SN 2006gy. The complex and unique 
spectral evolution of SN 2006gy will be discussed in a 
later paper, when more complete data are available. 

3.4. Initial Thoughts on Radioactive Decay and the 
Pair-Instability Hypothesis for SN 2006gy 

In previous sections, we have noted some obstacles, 
primarily observational in nature, with simple fireball or 
CSM interaction models as the engine for SN 2006gy. 
Although a suitable choice of extreme conditions may al- 
low them to work, at least in part, our observation of 
soft unabsorbed X-rays from SN 2006gy and the corre- 
sponding upper limits to the progenitor star's mass-loss 
rate make it worthwhile to consider other options. Pow- 
ering SN 2006gy with radioactive decay does not suffer 
from these problems, because this mechanism is known 
to work in other SNe. The question here centers around 
whether it is plausible to simply scale up the 56 Ni decay 
that powers fainter SNe, how that large mass of Ni may 
be created, and what happens to the radiation mecha- 
nisms in that extreme case. If SN 2006gy is powered 
by radioactive decay, the large Ni mass would require a 
pair-instability SN, as discussed below. 

Scannapieco et al. (2005) presented model light curves 
for pair-instability SNe, where the progenitor stars were 
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assumed to be red supergiants. The resulting light curves 
showed an initial small peak, but then a long, slow rise 
to maximum powered by 56 Ni and 56 Co decay. Some of 
their models get nearly the peak luminosity of SN 2006gy, 
but they rise more slowly to maximum than SN 2006gy 
did. However, their calculations were for zero metallic- 
ity, non-rotating stars with no pre-SN mass loss. Differ- 
ent assumptions about the metallicity, mass-loss, and the 
presence of rotational mixing may change things consid- 
erably (e.g., Maeder 1987; Yoon & Langer 2005; Woosley 
& Heger 2006). Also, if the progenitor of SN 2006gy had 
a small radius as we expect for an LBV (RSGs are not 
observed at high luminosity in normal- metallicity stars), 
then the initial peak may be lost due to adiabatic cool- 
ing, and the delayed rise after ^50 d would be dominated 
by 56 Co decay. Interestingly, this is similar to the case of 
SN 1987A, where the progenitor was a blue supergiant 
with a small radius, and where its late (70-100 d) peak 
was powered by radioactive 56 Co decay. SN 2006gy took 
a similarly long time to reach its peak luminosity, and 
its light curve thus far has a shape resembling that of 
SN 1987A (Fig. 2), but it was 250 times more luminous. 

In addition, the pair-instability models of Scannapieco 
et al. (2005) predict slow expansion speeds of ^5000 
km s _1 and the presence of H in the spectrum, again com- 
patible with SN 2006gy. These clues are tantalizing, and 
it would be interesting to see models for pair-instability 
SNe at metallicity closer to solar values and with com- 
pact progenitors. This is still somewhat virgin territory 
and will require continued observational constraints and 
detailed calculations to find a suitable model that will 
work for the case of SN 2006gy. Below we sketch a plau- 
sibility argument for the hypothesis that SN 2006gy was 
a pair-instability SN based simply on the required power 
source for its radiated luminosity. 

The i?-band magnitude at the peak of SN 2006gy was 
at least as bright as —21.8, but could have been signifi- 
cantly brighter because of our conservative assumptions 
for the reddening, as noted in §§2.1 and 2.2. Assuming 
no bolometric correction (again, conservative), this cor- 
responds to a peak luminosity of > (1.7 ± 0.3) x 10 44 
erg s . If this peak luminosity traces the instantaneous 
decay rate (Arnett 1982), we can estimate the necessary 
mass of initial nickel in the 56 Ni — > 56 Co — > 56 Fe de- 
cay. With a late peak at t ss 70 d, this will put us 
well into cobalt decay instead of nickel, as noted above. 
The radiated luminosity from cobalt decay (Sutherland 
& Wheeler 1982) is 

L = 1.42 x 10 43 erg s _1 e _t / ul d Af Ni /M Q (2) 

= 8 x 10 42 erg s" 1 M Ni /M , 

where Mni is the initial 56 Ni mass. The extreme luminos- 
ity of SN 2006gy, then, would require an extraordinarily 
high Ni mass of roughly 22 M Q to be synthesized in the 
explosion. This can be scaled down somewhat if CSM in- 
teraction contributes part of the energy, but unless that 
interaction dominates the light output, this large Ni mass 
cannot be explained with a core-collapse SN. (Compare 
this to a normal SN II arising from a star of 15-20 M Q , 
with a typical Ni mass of about 0.07 M Q .) 

The large Ni mass implicates a progenitor star that 
began its life with a mass well above 100 M Q . The con- 
sequences of this are potentially far-reaching, and could 



turn out to be the most interesting result of this study: 
namely, the only way to get such an extraordinarily high 
Ni mass to power the radiated energy would be from a 
pair-instability supernova, where the star's core is oblit- 
erated instead of collapsing to a black hole (Barkat et al. 
1967; Fraley 1968; Bond et al. 1984; Heger & Woosley 
2002). This type of supernova is only expected to occur 
in extremely massive stars. For the mechanism to work in 
the modern universe, even the most massive stars would 
need to retain most of their initial massive envelopes, 
providing a self-consistent interpretation of SN 2006gy in 
light of other evidence for its high mass discussed here. 
This is not wild speculation — it may even be the most 
promising explanation — but it deserves close scrutiny 
because of its far-reaching importance. 

As SN 2006gy continues to evolve, it will become eas- 
ier to determine if 56 Co decay or CSM interaction is 
the power source. If CSM interaction drives the visible 
light, we might expect the light curve to plummet precip- 
itously, down to the luminosity of a normal SN II, when 
the shock reaches the outer extent of the LBV shell. Such 
a drop occurred in SN 1994W, although the light-curve 
shape of SN 2006gy so far is quite different from that of 
SN 1994W (Fig. 2). On the other hand, if SN 2006gy con- 
tinues to decay smoothly from its peak, like SN 1987A 
but at an elevated luminosity, then it was almost cer- 
tainly a pair-instability SN event because of the large 
nickel mass required. So far, SN 2006gy shows no sign of 
plummeting — in fact, the latest photometry seems to 
imply that it is settling onto a plateau. 

Of course, SN 2006gy could be a combination of both 
CSM interaction and pair instability. Any very massive 
star capable of suffering a pair-instability SN is likely to 
have a strong stellar wind in its late pre-explosion stages 
anyway, consistent with the values of (1-5) xlO -4 that 
we infer from the X-ray interaction. The pair-instability 
SN models of Heger & Woosley (2002) predict mass-loss 
pulses that precede the final explosion. In fact, the ob- 
served optical spectrum of SN 2006gy requires that CSM 
interaction is occurring at some level, but the critical 
question is whether this interaction is capable of pow- 
ering the enormous continuum luminosity of SN 2006gy. 
Current indications are that it cannot. 

3.5. A Massive Circumstellar LBV Nebula 

Independent of the energy-budget arguments, the 
properties of the unshocked circumstellar gas around the 
progenitor of SN 2006gy are also consistent with the in- 
terpretation that it was a very massive star, and provide 
critical clues that strongly refute the hypothesis that it 
was powered by the Type la explosion of a low-mass star 
interacting with dense CSM. The high-resolution spec- 
trum in Figure [5] contains a narrow component to the 
Ha line, which also exhibits a clear P Cygni absorp- 
tion profile. It indicates that the SN is expanding into 
a hydrogen-rich dense stellar wind or outflowing circum- 
stellar nebula of the progenitor star, which has an expan- 
sion speed of 130-260 km s _1 indicated by the absorption 
component. This same narrow absorption component is 
seen in other lines in the spectrum of SN 2006gy, such as 
He I (He I A6680 and Fe II lines are shown in Fig. 4), Si II, 
Fe II, Ca II, O I, etc. The narrow He I lines are unusual, 
and may suggest He-enriched material in the CSM. 

This expansion speed is a critical clue to the nature 
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Fig. 7. — The long-slit Keck/DEIMOS spectrum of SN 2006gy and NGC 1260 in the region around Ha. It includes the central point 
source SN 2006gy at the zero-offset position, plus extended emission from the host galaxy NGC 1260 on either side of it. The extended Ha 
and [N n] emission, which follows the rotation curve of the galaxy and has an [N o]/Ha intensity ratio typical of H II regions, indicates 
that NGC 1260 does have active star formation. The light row below SN 2006gy is a bad row in the CCD and has been masked. 



of the progenitor star that cannot be neglected. It is 
much faster than typical wind speeds of AGB stars (10- 
20 km s _1 ), effectively ruling out the interpretation of 
SN 2006gy as a SN Ha. While it is unclear if the ex- 
pansion speed itself is in direct conflict with an inter- 
pretation involving a common envelope mass-loss phase 
(e.g., Taam & flicker 2006, and references therein), as 
suggested by Livio & Riess (2003) to explain the prop- 
erties of SN 2002ic, and in the first version of Ofek et 
al. (2007) to explain SN 2006gy, that interpretation is 
ruled out for SN 2006gy based on the energy budget (see 
§3.3). This speed is also too fast for a RSG wind (20-40 
km s _1 ), making it difficult to believe that the progen- 
itor star had an initial mass in the range 10-40 M Q . 
Moreover, the speed is an order of magnitude too slow 
for the wind of an O-type supergiant, H-rich WN, or 
Wolf-Rayet (WR) star progenitor. On the other hand, 
this speed is entirely consistent with an LBV wind or 
nebula (e.g., Smith et al. 2004; Smith 2006). Similar 
absorption speeds were seen in the narrow P Cygni ab- 
sorption of SN 1998S, which Fassia et al. (2001) also in- 
terpreted as a prior blue-supergiant phase. Chugai et al. 
(2002), however, interpreted it somewhat differently as a 
fast blue-supergiant wind sweeping into a red-supergiant 
wind. Significant acceleration of the slow red-supergiant 
wind would require a swept-up mass comparable to the 
fast-wind mass shortly before the SN, which makes this 
scenario implausible in the case of SN 2006gy because of 
the large mass implied. The typical LBV ejecta speed 
agrees well with our constraints from §3.3. 

Narrow blueshifted absorption components similar to 
Ha are seen in a number of other lines throughout the 
spectrum of SN 2006gy along with some relatively broad 
blueshifted absorption (Figs. [Hand [5]). Those absorption 
features are not always present (Fig. 3]) , while narrow Ha 
emission remains. Thus, we cannot be certain that the 
narrow emission and absorption components of Ha con- 
stitute a true P Cygni scattering profile, so we consider 
both cases here. For each case, the luminosity of the 
narrow Ha emission component is a relevant quantity. 
At a distance of 73 Mpc, the luminosity of the narrow 
emission component of the Ha line on day 96 (Fig. [5]) is 
Lhq ~ (1-3 ± 0.3) x 10 6 L (the absolute flux was cali- 
brated by scaling the red continuum to match observed 
KAIT photometry at the appropriate date and correct- 
ing for Ar ks 1 mag). Note that the true luminosity 
may be somewhat larger than this because the apparent 



luminosity may be reduced by the blueshifted narrow ab- 
sorption. 

If the narrow Ha component arises in an unshocked 
CSM wind, we can make a rough estimae of the den- 
sity immediately outside the radius of the shock, given 
by Rs = V s t, where we again take 14=4000 km s" 1 and 
t — 96 d is the time the Keck spectrum was taken. Such 
estimates are plagued with uncertainties in the ioniza- 
tion fraction and H mass fraction, so the estimate below 
is a lower limit assuming fully ionized pure H gas. Fol- 
lowing equation (1) of Chugai & Danziger (2003), for 
example, our measured value of Lhq implies a density 
of roughly 2xl0 8 cm" 3 just outside R s sa 3.3 x 10 15 cm 
(adopting a^ a = 8.64 x 10~ 14 cm 3 s _1 for the Case B 
Ha recombination coefficient as noted below). If V w is 
taken to be 200 km s -1 , this implies a mass-loss rate for 
the progenitor star of roughly 0.01-0.02 M Q yr" 1 . While 
this is an exceptionally high mass-loss rate, higher than 
what we infer from the X-ray emission (§3.2), it still falls 
short of what is required to power the visual luminosity of 
SN 2006gy by more than a factor of 10-20. It is interest- 
ing, however, that this value is comparable to progenitor 
mass-loss rates estimated for other SNe Iln with similar 
narrow Ha P Cygni features from the unshocked CSM, 
such as SNe 1997ab and 1997eg (Salamanca et al. 1998, 
2002). 

If the narrow Ha emission component arises instead 
from unshocked ionized gas in a detached CS shell neb- 
ula, however, then the mass implied would add yet an- 
other requirement that the progenitor star was very mas- 
sive. It may arise in a circumstellar shell like the Ho- 
munculus Nebula of 77 Carinae (Smith 2006), for exam- 
ple. Using Ljjai and assuming that the line originates 
from a circumstellar shell nebula of constant density, the 
ionized gas mass can be expressed as 



Mi 



Ha 



hva 



eff 



c rr 



(3) 



8.64 x 



where hv is the energy of an Ha photon, a' 
10~ 14 cm 3 s _1 is the Case B Ha recombination coeffi- 
cient, and n e is the average electron density. This yields 
Miia ~ 11-4 M (Ln a /n e ). We do not know the electron 
density in the nebula around SN 2006gy, but values of 
10 5 -10 6 cm -3 are the highest densities typically seen in 
young LBV nebulae like the one around 77 Carinae (Smith 
2006). With the observed Ha luminosity and densities 
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of this order, the nebular mass is probably above 5 M@, 
and it could plausibly be as high as 20-30 Mq. Lower 
densities typically seen in circumstellar nebulae around 
lower-mass stars would require implausibly high emitting 
masses to account for the observed radiation, exceeding 
their own stellar masses. Environments this massive ob- 
viously cannot be produced by low-mass stars and are not 
seen around moderately massive stars of 20-40 M , but 
they are quite typical of the nebular shells around LBVs 
with L > 10 6 L (Smith & Owocki 2006), which descend 
from stars with initial masses of 80-150 Mq. Such large 
masses are consistent with the >12.5 M© nebula around 
r] Car (Smith et al. 2003). 

Thus, the flux of the narrow Ha component that we 
observe is only likely to arise in the circumstellar nebula 
of an extremely massive star. Taken together, this high 
mass and the shell's expansion speed give self-consistent 
evidence that the progenitor star was indeed very mas- 
sive. This line of reasoning is independent of the uncer- 
tainty associated with the mechanism that powers the 
radiated energy of the SN. It is also consistent with the 
presence of strong hydrogen lines in the spectrum, since 
LBVs have not yet shed their H envelopes. Although 
dominated by hydrogen, LBV shells also have elevated 
helium abundances, consistent with the presence of nar- 
row He I lines in the spectrum of SN 2006gy. If SN 2006gy 
really is surrounded by a dense LBV nebula like that of 
•q Carinae, then we might expect to see strong, narrow 
emission lines of [N n] AA6548, 6583 in its late-time spec- 
tral evolution, since LBV nebulae like that of r\ Car tend 
to be enriched with CNO-cycle ashes (Smith & Morse 
2004). 

3.6. Do We Expect Massive Stars in the Host Galaxy? 

SN 2006gy has been compared (Ofek et al.) to two 
peculiar supernovae, SN 2002ic and SN 2005gj (Fig. 2), 
which have been proposed as SNe la interacting with 
dense CSM (the so-called "Type Ha" SNe) as noted ear- 
lier. One factor that motivated Ofek et al. (2007) to 
originally favor the SN Ha hypothesis was that the host 
galaxy, NGC 1260, was apparently not a star-forming 
galaxy. It should not have massive stars, because SO 
galaxies are dominated by old stellar populations. 

We note, however, that the SN host is actually a pecu- 
liar SO/Sa galaxy with infrared (IR) emission from dust. 
NGC 1260 was detected by the Infrared Astronomical 
Satellite (IRAS), and Meusinger, Bruzendorf, & Krieg 
(2000) give an infrared luminosity of log(Lj#/L©) = 
9.85. According to Kennicutt (1998), this would trans- 
late to a star-formation rate of ~1.2 M© yr^ 1 , which 
is certainly high enough to permit this galaxy to host 
some massive young stars. Furthermore, we detect ex- 
tended Ha and [N ii] A6583 emission from the galaxy 
in our spectra; Figure [7] presents the original long-slit 
Keck spectrum before the Ha profile of SN 2006gy was 
extracted, revealing extended emission from gas that fol- 
lows the rotation curve of the host galaxy. These emis- 
sion lines, having intensity ratios typical of H II regions, 
are indicative of current star formation and are absent in 
non-star- forming galaxies. 

A related point concerns the statistics involved. 
SN 2006gy is the most luminous SN seen to date, but 
it is also spectrally peculiar, almost in a class by itself. 
Its unusual nature would not be at all surprising, in prin- 



ciple, if its origin were the explosion of a >100 M© star, 
since these stars are so phenomenally rare to begin with. 
On the other hand, if it results from normal evolution 
for low-mass stars or even moderately massive stars of 
10-40 M , then we would expect such events to be more 
common. 

4. SUMMARY: EXPLOSION AS A MASSIVE LBV AND THE 
RELEVANCE OF A PAIR-INSTABILITY SUPERNOVA 

All available observations are broadly consistent with 
the hypothesis that the progenitor of SN 2006gy was a 
very massive star that retained a massive hydrogen enve- 
lope until it exploded. Retaining this envelope does not 
mean that the progenitor was a RSG; the most luminous 
stars evolve to the LBV phase before losing their en- 
velopes, and during that phase they are hot supergiants 
with relatively small radii. This can strongly affect the 
early light-curve shape. A mass below 60 M may be 
possible if the event was powered by CSM interaction, 
but then one must invoke exceptional conditions incon- 
sistent with observed properties of stars below that mass. 
If CSM interaction dominates, we find it more likely that 
the progenitor star had an initial mass of 100-150 M , 
although we still lack a satisfactory explanation for the 
weak unabsorbed X-rays in that case. 

By contrast, the huge radiated luminosity, the long du- 
ration, the presence of hydrogen in the spectrum, the low 
expansion speed of the SN ejecta, and the various critical 
clues from the circumstellar environment are all consis- 
tent with the hypothesis that this event was powered by 
a pair-instability supernova that also has some moderate 
CSM interaction, implying that the progenitor star's ini- 
tial mass may have been near the upper mass limit for 
stars of ^150 M (Figer 2005). Regardless of the power 
source, several clues hint that the progenitor star may 
have resembled the LBV star r\ Carinae. 

If this hypothesis of explosion as a massive LBV is cor- 
rect, it would have important consequences for our un- 
derstanding of stellar evolution. It is currently thought 
that variability in the LBV phase is responsible for the 
mass shedding that marks the transition from the end of 
core H burning to core He burning, after which a star 
appears as a He-rich WR star (Abbot & Conti 1987; 
Langer et al. 1994; Smith & Owocki 2006; Smith et al. 
2004). During this brief evolutionary phase, a massive 
star might undergo sequential bursts of mass loss when 
it can repeatedly shed more than 10 M of material in a 
decade (Smith & Owocki 2006). These events are seen in 
other galaxies as faint SNe Iln, or "supernova impostors" 
(Van Dyk 2005, and references therein). They may dom- 
inate the mass loss of the most massive stars, shedding 
more total mass than line-driven winds during the star's 
lifetime (Smith & Owocki 2006). Consequently, LBV 
stars are frequently surrounded by circumstellar nebu- 
lae with masses of order 10 M©, like the one that may 
reside around SN 2006gy. It would appear that one of 
these events may have occurred within a decade or so 
immediately preceding SN 2006gy. 

The core He burning WR phase that should follow af- 
ter the massive hydrogen envelope is stripped away is 
expected to last a few hundred thousand years before 
the star reaches even more advanced stages of nuclear 
burning and finally explodes (Abbott & Conti 1987). If 
LBVs explode before reaching the WR phase, though, 
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it means that they could be in more advanced stages of 
nuclear burning than currently predicted by stellar evolu- 
tion theory. SN 2006gy adds to mounting evidence (e.g., 
Smith & Owocki 2006; Kotak & Vink 2006; Gal- Yam 
et al. 2007; Smith 2007) that stars may explode "early" 
during the LBV phase, and it hints that reaching the pair 
instability could be a reason for this. 

It seems intuitively possible, although difficult to 
prove, that it would be the most massive LBVs above 
~100 M Q that are more likely to explode prematurely 
as they have a greater burden of removing their massive 
envelopes before transitioning to WR stars. Gal- Yam 
et al. (2007) have drawn a similar conjecture, consider- 
ing LBVs as the most likely progenitors of SNe Iln. If 
the most massive stars can indeed explode before the 
WR phase, then our current ignorance of the instability 
underlying the LBV phase presents a critical challenge. 
The possibility that SN 2006gy could have been a pair- 
instability supernova weighs heavily upon the impor- 
tance of understanding these LBVs as well. SN 2006gy 
may be giving us a clue that the wild instability of the 
most luminous LBVs like r\ Carinae could be early warn- 
ing signs of a massive star's imminent demise, and there 
may be theoretical reasons to think this is the case. One 
implication is that we had better keep a watchful eye on 
r\ Carinae. 

The chief reason why pair-instability SNe are expected 
to occur for high-mass stars in the early universe is be- 
cause their low metal content is expected to reduce their 
mass-loss rates, causing them to retain their massive H 
envelopes (Heger et al. 2003; Heger & Woosley 2002; al- 
though see Smith & Owocki 2006). Also, the initial mass 
function of the first stars is thought to have been skewed 
to higher masses due to the lack of metal cooling and 
consequent fragmentation in the star-formation process 
(e.g., Bromm & Larson 2004). SN 2006gy may have been 
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a very massive star that exploded as an LBV before it 
could shed its H envelope, and it may have done so by 
the pair-instability mechanism. 

The fact that SN 2006gy was able to explode success- 
fully instead of winking away into a black hole has far- 
reaching implications. In particular, one primary goal 
of the James Webb Space Telescope will be to search for 
these first explosions in the universe, and the brilliant 
display of SN 2006gy may bode well for the possibility of 
their infrared detection at high rcdshift. 
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